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News on Disconnected Diagrams* 
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b Physics Department, University of Wuppertal, D-42097 Wuppertal, Germany. 

We present evidence for disconnected contributions to the ovjv term and the flavor singlet axial coupling of 
the proton on full QCD configurations, which are obtained by means of improved stochastic estimator techniques. 

Furthermore we discuss results from the fermionic determination of the topological charge of the QCD vacuum 
(in the spirit of the Athiya-Singer theorem) again achieved with stochastic estimator methods. It turns out 
that this approach provides an monitor for the tunneling efficiency of the HMC on QCD with dynamical Wilson 
fermions which is independent on the pure gluonic method that implies cooling. 
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1. Introduction 

Flavor singlet matrix elements, like the ovat 
term or the axial coupling of the proton g A , have 
been widely discussed in the past. Recently first 
lattice calculations worked out the connected as 
well as the disconnected diagrams contributing to 
oVat and g° A 0J3]. It turned out that the discon- 
nected diagrams define a severe computational 
problem whereas the methods for the connected 
amplitudes work reasonable good. Two differ- 
ent approaches to tackle the pure vacuum fluctu- 
ations of disconnected diagrams were suggested 
in the literature, namely the wall-source method 
|3| and stochastic estimator techniques (SET) [Q. 
In the following we focus on SET with Z2 noise, 
coming back to the wall-source later on. 

2. Stochastic Estimator Techniques 

2.1. Standard Z2 Noise 

Stochastic matrix inversions are based on a set 
of appropriate source vectors r\ where all compo- 
nents carry a Z2 noise and fulfill 



(m) = o, 



(ViVj) 



(1) 



in the stochastic average. For infinite number of 
estimates M _1 can then be obtained from 



(2) 



where x solves Mx = rj for different source vectors 
r\. The number of estimates equals the number of 
matrix inversions and the accuracy is bounded 
according to the computer facilities available. 

2.2. Improving SET 

The stochastic errors for off-diagonal traces, i.e. 
Tr^M -1 ), Tr (71 75 A/ increase drastically 
and improved SET methods become inevitable. 
Hence, we worked out the effective noise reduc- 
tion of exact inversion in the spin subspace of 
Wilson's fermion matrix M |^J. However, this 
method reduces the number of estimates for a 
constant number of inversions by a factor 4, and 
could influence the accuracy of Tr(M~ 1 ) and 
related observable. To cope with this problem 
diagonal-improved SET was invented. 

The error of M~ x for finite numbers of esti- 
mates, 



(3) 



collects a large contribution when k = i due to 
M~l X >> for the Wilson fermion matrix. 

The error A can be suppressed by subtracting the 
current estimate for M^ 1 after each inversion: 

E[M^\ = M^ + ^M^inw) 



(4) 
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where i ^ j. This diagonal-improved SET reduces 
the standard error for Tr^M -1 ), Tr (71 7sM~ 1 ) 
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and Tr(727 5 A/- 1 ) about 40% without addi- 
tional matrix inversions and does not reflect on 
Tr^M- 1 ). 

3. Flavor Singlet Matrix Elements 

We look at matrix elements of the proton 
(p\qTq\p) = Aq at zero momentum when qTq is 
a flavor singlet bi-quark operator, i.e. T = 1 for 
the a^pf term and r = 75 7 M for the axial coupling 
g° A . The ratio of the three-point function and the 
proton propagator 



R(t) 



(pW(SngTg)P(O)) 

<P(*)P(0)> 



const + Z r 1 (p\qTq\p) t , 



(5) 



is for large t proportional to Aq ||. Connected 
diagrams are calculated with the conventional 
source method whereas for the disconnected 
amplitudes we applied improved SET to estimate 
Tr(rM _1 ) which is required in eq. §. The 
SESAM gauge-field configurations are generated 
with 2 flavors of dynamical fermions on 16 3 x 32 
lattices with the standard Wilson action (j3 = 5.6, 
K = 0.1560-0.1575). Figure |l| shows R(t) on 200 
configurations at k = 0.156 for ovjv and g A , re- 
spectively. For scalar insertions the disconnected 
amplitudes are compatible to the connected ones. 
Au coni Ad con and Aq dls = Au,d dls are positive. 
Axial connected insertions have opposite sign and 
Aqdis is rather small and negative. Even with im- 
proved SET the errors for Aqdis are still around 
30% . The outstanding calculation of As^is for 
the strange quark will complete the analysis for 
a n N and g A and is in advanced progress. Finally 
we should point out that the wall-source method 
is too noisy for the axial disconnected insertion 
and improved SET allows the determination of 
g A on our set of full QCD configurations. 

4. Topology and SET 

Topological properties of the QCD vacuum are 
essential for the understanding of flavor singlet 
physics, e.g. the unexpectedly large 77' mass or 
the proton spin. Here the topological charge Q 
plays a crucial role since the Witten-Veneziano 




Figure 1. R(t) and R 9A (t) for the scalar 
and axial insertion on 200 configurations, K = 
0.156. The linear slopes are proportional to 
Au con , Ad con and Aq, 
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formula relates the topological susceptibility, x = 
(Q)/V, directly to the rj' mass. Instead of us- 
ing the more usual field-theoretic H or geometric 
definitions || of Q, we look at the Athiya-Singer 
index theorem in the continuum 



Q 



m q Tr('j 5 G(x, x')) 



(G) 



where are the numbers of left and right 



m q the quark mass 
On the lattice this 



handed chiral eigen-modes, 
of the propagator G(x,x'). 
turns into 

Ql = m^pTr^hr 1 ). (7) 

Here M ~ 1 is the inverse fermion matrix and 
renormalization constant fllpf . Again we used im- 
proved SET to estimate the right hand side of eq. 
0. In figure || the results for on a set of HMC 
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trajectories from SESAM lattices (k = 0.1575) 
are displayed. The topological charge via the in- 
dex theorem is compared to Ql when we use cool- 
ing and the field-theoretical definition [||. The 
best overlap of both data is achieved by setting 
rriqKp = 2m q in eq. [7] and a further determina- 
tion of Kp is imperative. The strong correlation 
of the data indicates that the index theorem re- 
mains valid on the lattice. Finally two completely 
different methods display the tunneling efficiency 
of the HMC algorithm through various topologi- 
cal sectors and justify the ongoing calculation of 
X, the rj mass and the proton spin. 
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Figure 2. The field-theoretical definition of the 
topological charge Ql after cooling and via the 
index theorem on a set of HMC trajectories for 
k= 0.1575. 



5. Conclusion and Outlook 

We presented an improved stochastic estimator 
technique and found evidence that disconnected 
diagrams contribute to the ovjv term and to the 
flavor singlet axial coupling of the proton, g A . 
It turns out that the disconnected contribution 
Ait, ddis to ovjv is comparable to the connected 
amplitude, whereas Ag^is for g A is rather small 
but non- vanishing. The final analysis for both 
flavor singlet operators, o^n and g A , will be pub- 
lished soon. 



Furthermore we achieved values for the topo- 
logical charge Ql with stochastic matrix inver- 
sion. These data are in good agreement with the 
results for Q l when the field-theoretical definition 
and cooling is applied. Both methods demon- 
strate the tunneling efficiency of the HMC on 
SESAM lattices throughout different topological 
sectors. 
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